Manifestation of chiral tunneling at a tilted graphene pn junction 
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Electrons in graphene follow unconventional trajectories at PN junctions, driven by their pseu- 
dospintronic degree of freedom. Significant is the prominent angular dependence of transmission, 
capturing the chiral nature of the electrons and culminating in unit transmission at normal incidence 
(Klein tunneling). We theoretically show that such chiral tunneling can be directly observed from 
the junction resistance of a tilted interface probed with separate split gates. The junction resistance 
is shown to increase with tilt in agreement with recent experimental evidence. The tilt dependence 
arises because of the misalignment between modal density and the anisotropic transmission lobe 
oriented perpendicular to the tilt. A critical determinant is the presence of edge scattering events 
that can completely reverse the angle-dependence. The absence of such reversals in the experiments 
indicates that these edge effects are not overwhelmingly deleterious, making the premise of transport 
governed by electron 'optics' in graphene an exciting possibility. 

PACS numbers: 72.80.Vp, 73.63.-b, 72.10.-d 



A striking feature of electron flow in graphene, gated 
uniformly or electrostatically 'doped' into junctions, is 
the non-trivial dynamics of its pseudospins arising from 
its orthogonal dimer basis sets. The overall photon like 
dispersion propels electrons along trajectories intuitive 
of Snell's law, conserving quasi-momentum components 
transverse to the interface. However, the corresponding 
'Fresnel equations' are qualitatively different from their 
optical counterpart, determined by conservation of pseu- 
dospins. In particular, graphene electrons are chiral in 
nature meaning the pseudospin components are related 
to the direction of momentum. This results in perfect 
transmission at normal incidence [l] regardless of voltage 
gradients across the junction (Klein tunneling). For other 
incident angles, the spinor mismatch leads to a unique 
angle dependent transmission across the junction. Thus, 
while conventional electronics in graphene faces possibly 
steep challenges [2l, the dynamics of pseudospintronics 
can usher in novel concepts such as electronic Veselago 
lens [3 , subthermal switches driven by geometry induced 
metal-insulator transition [4 and Andreev reflections [5 . 

Despite the exciting physics of chiral electron flow, its 
measurable signatures have so far been sparse and indi- 
rect. Signatures of Klein tunneling were seen [6l [7] in 
the preferential transmission of normally incident carri- 
ers predicted in [8 . A more direct measurement was 
the conductance oscillation in an npn structure [9 . The 
reflection amplitude undergoes a phase shift of tt at nor- 
mal incidence under the action of a magnetic field, due to 
the cyclotron bending of the carriers [10 . However, the 
main underlying physics of the angle dependent electron 
transmission has not been explicitly measured. Neither 
has there been a proper model that can capture both 
the quantum mismatch of spinors over the entire dop- 
ing regime, as well as diffusive scattering to explore their 
robustness to impurity and edge scattering events. 

In this paper we focus on a tilted GPNJ (i) to show that 
it serves as an explicit signature of chiral tunneling. We 
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FIG. 1. (Color online) (a) false color SEM image 11 and (b) 
device schematic. Two split gates, separated by a distance d 
create pn junction with linear potential variation across the 
junction. The gates are placed at an angle 5 with respect to 
the transport direction, (c) Experimental [11 , (d) theoretical 
junction resistance [Rj) vs volt age /energy for 5 — 0,45°. In 

(d), we plot Rj vs Ef — Tivf\/ TvacCc \VGi\/q from atomistic 
NEGF calculation for a lOOnm wide graphene sheet. 

show that the junction resistance (similar to the odd re- 
sistance shown in other experiments ^l7j) is higher than 
the non-tilted device (Fig. 1). We argue that this en- 
hancement originates from the chiral nature of graphene 
electrons which manifests itself through the highly angle 
dependent transmission characteristics of GPNJ (Fig. 2). 
The angular transmission lobe, oriented perpendicular to 
the interface is rotated with the tilt, where fewer trans- 
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mitting modes exist. Therefore the conductance modu- 
lation would not happen for non-chiral, non-relativistic 
electrons with isotropic transmission. The results follow 
closely with our recent transport measurements of a tilted 
GPNJ in a structure that has separately controlled split- 
gate voltages [m [12] . (ii) We present an exact analyt- 
ical solution to the spinor mismatch problem as well as 
non-equilibrium Green's function (NEGF) based atom- 
istic numerical calculation. An efficient matrix inversion 
algorithm is employed to reach near experimental dimen- 
sions and capture both quantum mechanical and diffu- 
sive contributions to the overall resistance. We find that 
charged impurity scattering dilutes, but does not elim- 
inate the modulation in conductance (Figs. 1, 3). (iii) 
Notably, we demonstrate that multiple scattering events 
at edges can reverse the trend in modulation, giving an 
interfacial resistance that decreases with tilt (Fig. 4). 
Such decrease has been seen in the past [13 , but their 
physical origin has not been identified so far. The absence 
of such a reversal in experiments surprisingly points to 
their elimination, possibly through incoherent scattering 
processes dominant at the strained edges. 
Analytical model. The conductance of a GPNJ can be 
written as, G = Go'^f^ T{Oi) where Go = ^q^/h is the 
conductance quantum including spin and valley degen- 
eracies and T{Oi) is the incident angle dependent trans- 
mission probability with Oi = tan~^{ky/kx). M is the 
number of modes from the incident side for a given Fermi 
energy Ep relative to its Dirac point and can be approx- 
imated as M WIEfI/ttTivf over the linear E-k regime. 
The angle dependent transmission T is obtained by pseu- 
dospin conservation across the junction |4] 
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This is a general form of the transmission expression in 
[HI [14] and works for the entire voltage range from PN to 
A^A^ junction. The incident and refracted angles Oi and Or 
are related by Snell's law [3 , Ep sin^^ = {Ef — Vq) sin^^, 
with Vq (Fig. 1) being the voltage barrier across the 
junction, and kp = Ep/fiVF' The Snell's law arises from 
transverse quasi-momentum conservation (energy band 
diagram in Fig. 1(b)). For the rest of the paper, we use 
the average transmission of all modes, Tav = G/M. 

For a tilted junction, the incoming mode angles are 
modified, so that the conductance becomes, 



M{Ef) 

G{Ef)=Go T{Ef,0, 



(2) 



where S is the tilt angle as shown in Fig. 1. In Eq. 
[2] the effective split between the two gates is d/cosS. 
As a result of the angle dependence, the transmission 
lobe at a particular energy will rotate by the tilt angle 
(Fig. |2|l). The transverse wavevector ky = ksinO gives 
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FIG. 2. (Color online) (a) Angular transmission for various 
tilt angles, (b) With tilt, the transmission lobe moves into a 
low mode density (^ cosO) area giving (c) a gradual decrease 
in transmission (Eq. [3| for a symmetric pn junction, (d) junc- 
tion resistance at OK predicted from Eq. |4] We see resonances 
which become more pronounced as we go into smaller systems. 



A6> = ISkyjKcosO^ so that the mode density decreases 
as we go to higher angles relative to the transport axis. 
A tilt at the junction thus shifts the transmission win- 
dow onto a high angle region where the mode density is 
less, decreasing the overall transmission (Fig. In the 
limit when the no. of modes is very few, the experiment 
will give the mode resolved angle dependent transmis- 
sion properties (Fig. 2(a)). For an abrupt, symmetric pn 
junction, the transmission expression reduces to cos'^0 
from Eq. [l]and it is easy to see the impact of tilt, 
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The factor 2/3 arises from the wavefunction mismatch 
across the junction and the tilt introduces an additional 
scaling factor, which further reduces conductance. The 
gradual decrease of Tav with S in Fig. |2|c) constitutes a 
direct manifestation of chiral tunneling in graphene. 

To connect with experimental measurements, we next 
analyze the variation of the junction resistance in pres- 
ence of an intrinsic background doping (Vdp) in the 
graphene sheet (Fig. 2d). We vary the gate voltages so 
that Vgi = —Vg2 but a nonzero Vdp makes it an asym- 
metric GPNJ. The effective gate voltages on the graphene 
sheet are 0^0(^01 -\-Vdp) and 0^0(^02 + ^dp), where ac 
is the capacitive gate transfer factor. The junction resis- 
tance can be written as [15 , 



(4) 
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Fig.[2[d) plots Rj against Ef = Tivp \/ (T^)(y^GCG \VGi\/q, 
the amount of shift in the Dirac point by Vqi-, for a 
lOOnm wide graphene sheet with a spht gate separation 
d = 200nm. For the voltage range \Vgi\ < Vdp, we are 
in the n+n regime for positive Vdp and p'^p for neg- 
ative. Under these near homogeneous conditions, the 
junction resistance predicted by Eqs. 1-4 is small, be- 
cause the pseudospin states match and there is no WKB 
tunneling term in Eq. [l] When \Vgi\ > Vdpi we are in 
the pn junction regime and resistance jumps to a high 
value, primarily due to the WKB factor in Eq. [l] (sim- 
ilar trend was seen in 00). The rate of change in Rj 
with Vqi is determined by d. For the 45^ tilted junction, 
the junction resistance is higher than the non-tilted re- 
sistance. We see oscillation in the resistance (Fig. 2(d)) 
for the single np junction. This is different from the in- 
terference oscillation in Ref. [9 for the resonant cavity 
formed in npn strucure. This can be understood from 
the conductance in the pn junction regime, simplified as 
G{kF) ^ GoYli^ ^^'P{~'^^Fdsin'^0i/2) large d. With 
increasing gate voltage (higher /c^?) we have more modes 
(M) in the summation with each mode transmitting with 
exponentially reduced magnitude. The two opposing ef- 
fects generate a sequence of peaks and valleys and dom- 
inate when we sum over a few modes (either with quan- 
tization or tilt) at low temperature. 

Numerical model. An atomistic, numerical calculation 
of the junction resistance is shown in Fig. 1 (d) at 80K 
temperature. We use NEGF formalism for a lOOnm wide 
graphene sheet with d = 200nm, close to experimental di- 
mensions (width ~200-300nm). A single Pz orbital basis 
for each carbon atom is used to compute the Hamiltonian 
while the contact self-energies Si ^2 are calculated us- 
ing a recursive technique. The retarded Green's functions 
are calculated as = {EpI - + V - Si - T.2)~^ us- 
ing the algorithm in [16 , V is the electrostatic potential 
inside the device. In units of 2q^ /h^ the conductance is 
calculated as G = riG^r2G^, where the contact broad- 
ening functions F are the anti-Hermitian components of 
S. For ballistic transport, G equals M, the number of 
modes for a uniformly gated sheet, and MTav for a pn 
junction. Combining these two, we calculate Tav and 
Rj from Eq. |4j which shows a jump with tilt in the pn 
junction regime, very similar to the experiment. 

Charged impurity scattering dilutes tilt dependence. 
The experimental device is on a Si02 substrate and the 
transport is diffusive with a mobility varying from 700- 
3000 cm^/V-s. It is natural to inquire how the the- 
oretical model, which so far had no scattering, corre- 
sponds to experiments. To explore this feature, we in- 
cluded the impact of charged impurity scattering in our 
model. We use a sequence of screened Gaussian po- 
tential profiles for the impurity scattering centers [17]- 
U{r) = EnLTUnexp{-\r-rn\y2e), that speci- 
fies the strength of the impurity potential at atomic site 
r, with rn being the positions of the impurity atoms and 
C the screening length (^ 8 times the C-C bond for long 
range scatterers). The amplitudes Un lie in the range 
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FIG. 3. (Color online) Impact of charged impurity scattering, 
(a) conductance asymmetry is diluted due to impurity poten- 
tials, the ballistic resistance is normalized for comparison, (b) 
Reduced asymmetry results in lower junction resistance for 
both tilted, non-tilted devices, thus retaining their difference. 



[—(5, 6] (^ 0.5 times the C-C coupling parameter) and 
Nimp is the impurity concentration (~ 5 x 10^^ cm~^). 
Note that the purely diffusive model discussed in [6|[7[[2Q] 
ignores the quantum mechanical spinor mismatch and 
WKB scaling and therefore underestimates the junction 
resistance for cleaner samples. Our NEGF based numer- 
ical model on the other hand captures both the quantum 
mechanical and impurity limited resistance contributions 
simulataneously. The junction resistance is now calcu- 
lated by eliminating the contact and device resistance 

m 

Rj = [R{Vgi,Vg2) + R{Vg2, Vgi) 

-R{Vgi,Vgi) - R{Vg2, Vg2]/2, (5) 

where the first two terms contain the junction resistance, 
while the last two don't. 

Fig. 3(a) shows the impact of the impurity scatter- 
ings on the total resistance and 3(b) on the junction 
resistance. We take the average resistance over many 
impurity configurations. This puts a constraint on the 
computation size, so we show calculations this time for a 
smaller device (50nm wide). We find that both tilted and 
non-tilted junction resistances are suppressed, thereby 
retaining the difference between the two. This reduction 
in junction resistance with scattering is quite consistent 
with the experiment (Fig. 1(c), red line is for a 45^ de- 
vice with mobility 2270cm^/V-s, while the blue triangle 
is for 45^ with lower mobility, 700cm^/V-s). 

The reduction in the junction resistance from ballistic 
to diffusive transport can be understood from the trend 
in total resistance, shown in Fig. [sj^a). Now we keep 
Vgi fixed and vary Vg2 so that we go from n+n to np 
junction. We see a clear asymmetry in the R-Vg [3 fT4) 
(purple line normalized to the orange line in Fig. 3(a) 
for comparison). The asymmetry confirms the presence 
of pn junction, which reduces the conductance due to 
spinor mismatch. The presence of impurity scattering 
reduces this asymmetry while increasing the overall re- 
sistance (red line). The impurity potentials create a ran- 
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FIG. 4. (Color online) (a) Increase in conductance in a tilted 
GPNJ due to Edge Scattering (ES) in contrast with Fig. 2(c). 

(b) corresponding decrease in junction resistance due to tilt 

(c) Mechanism of edge enhanced conductance for a tilted junc- 
tion from atomistic NEGF calculation: some reflected elec- 
trons come back at the junction after edge reflection. 



dom potential variation throughout the graphene sheet 
on top of the applied gate voltages, thus blurring the 
presence of a pn junction. Therefore the resistance due 
to spinor mismatch becomes less noticeable (Fig. 3(b)). 
Indeed, the experimental data of the total resistance in- 
dicates an increase in asymmetry in the tilted junction 
J2], signifying an increase in the junction resistance. 

Edge scattering can reverse tilt dependence. A strik- 
ing feature on the experimental results is their agree- 
ment with Eq. |3] This match is remarkable, consid- 
ering that the equation was derived assuming no edge 
reflections and the fact that past numerical study |13j 
showed in fact an increase in conductance with tilt. We 
argue that the above reversal of junction conductance 
with tilt is entirely due to edge scattering events. In- 
deed, from an atomistic NEGF calculation with shorter 
widths than lengths, we find that the transmission now 
shows a pronounced local maximum (Fig. |4ja) orange 
line) in agreement with [13 , increasing thereby the junc- 
tion resistance. We summarize this in Fig.[4]^b), where an 
increasing tilt makes the resistance increase for the short 
channel 125nm x 50nm device (a transition from pur- 



ple circle to black diamond), but decrease for the longer 
200nm x 50nm device (orange square to red triangle). 
Bearing in mind that the gate split is lOOnm, the short 
channel device significantly reduces edge scatterings. 

To better understand the origin of such a resistance re- 
versal, we inject electrons with a small contact at the left 
edge (bright red spot in Fig. ^) and plot the spatial cur- 
rent density under a small drain bias. The numerically 
computed electron trajectories show how a tilt can en- 
hance forward scattering events at the edge and thus an 
increase in conductance. The enhancement arises from 
simple 'geometrical optics' dictated by Snell's law. We 
can identify the incident wide angle modes {6 > ir/A — S), 
for which the reflected 'ray' hits the upper edge with a 
positive X directed velocity. Such a mode will reflect back 
towards the junction again. The contribution from the 
positively directed edge scattering event is given by 
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Note that only the incident angles below the critical angle 
are considered while setting the limits of the integration. 

With the added contribution from edge scattering, the 
net mode- averaged transmission T = G/M is given by 



^total ^ T + (1 - T)Tedger] 



(7) 



where the Ts are extracted from the corresponding G/M 
ratios and 77 is a parameter that describes the efficiency of 
edge scattering. In the absence of edge scattering (77 = 0), 
Ttotai = T and decreases with tilt (Fig. 2c). However 
in the presence of strong edge scattering {r] = 1), the 
added forward edge scattering term in Eq. [7| closely re- 
produces the NEGF result with the local transmission 
maximum (Fig. |4|i, black dotted line) . Comparing these 
results with experiments indicates that such edge scat- 
tering events are clearly minor. We conjecture that the 
coherent forward scattering processes captured by NEGF 
can be diluted down in the experiments by the presence 
of incoherent scattering processes arising at the strained 
and rough edges of the graphene samples that tend to 
dephase or perhaps even trap the electrons. 

The experimental observation of chiral tunneling, par- 
ticularly in the face of impurity and edge scattering, 
opens up the possibility of graphene's 'geometric optics' 
based applications, e.g. lens, switches and wave-guides. 
Analogous results are expected in bilayer graphene, but 
not in achiral materials like 2-D hexagonal boron nitride. 
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